Evidence for nodeless superconducting gap in NaFei^Co^As from low-temperature 

thermal conductivity measurements 
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The thermal conductivity ol optimally doped NaFeo.972Coo.028 As (T c ~ 20 K) and overdoped 
NaFeo.925Coo.075 As (T c ~ 11 K) single crystals were measured down to 50 mK. No residual linear 
term kq/T is found in zero magnetic field for both compounds, which is an evidence for nodeless 
superconducting gap. Applying field up to H — 9 T (« H C 2/4) does not noticeably increase ko/T 
in NaFei.972Coo.o28As, which is consistent with multiple isotropic gaps with similar magnitudes. 
The ko/T of overdoped NaFe1.925Coo.075 As shows a relatively faster field dependence, indicating 
the increase of the ratio between the magnitudes of different gaps, or the enhancement of gap 
anisotropy upon increasing doping. 



PACS numbers: 74.70.Xa, 74.25.fc 

The discovery of high-temperature superconductivity 
in iron pnictides 1,2 has attracted great attentions. Many 
efforts have been devoted to determine the symmetry 
and structure of their superconducting gap,— which is 
one of the keys to understand the electronic pairing 
mechanism^ 

For the most studied (Ba,Sr,Ca,Eu)Fe 2 As 2 (122) 
system, while multiple nodeless superconducting 
gaps were demonstrated near optimal doping, 5-9 
nodal superconductivity was found in the extremely 
hole-doped KFe2As2 , 1Ql11 and highly anisotropic 
gap^ or isotropic gaps with significantly different 
magnitude o 12 ' 13 were suggested in the heavily electron- 
doped Ba(Fei_ 2 ,Co a; )2As2. Intriguingly, nodal su- 
perconductivity was also found in isovalently doped 
BaFe 2 (As 1 _ x P :r )2 and Ba(Bfei_ B Ru a .) 2 As2^il 

For the two stoichiometric LiFeAs and LiFeP (111) 
compounds, it was clearly shown that the former 
has nodeless superconducting gaps ) 18-21 and the lat- 
ter has nodal gap. 22 However, for another 111 com- 
pound NaFei-^Co^As, recently there were controver- 
sial experimental results on its gap structure . 23 ! 24 High- 
resolution angle-resolved photoemission spectroscopy 
(ARPES) measurements on NaFeo.95Coo.05As (T c ~ 18 
K) unambiguously showed nearly isotropic supercon- 
ducting gaps on all three Fermi surfaces,— but London 
penetration depth measurements claimed nodal gap in 
NaFei-zCozAs from underdoped to overdoped range (x 
ranging from 0.02 to 0.10).— To resolve this important 
issue, more experiments are highly desired. 

The ultra-low-temperature heat transport measure- 
ment is a bulk technique to probe the gap structure 
of superconductors.— Whether there is a finite residual 
linear term kq/T in zero field is a good judgement on 
the existence of gap nodes or not. The field dependence 
of kq/T can further give information of nodal gap, gap 
anisotropy, or multiple gaps. 25 During past three years, 



thermal conductivity measurements have provided valu- 
able experimental results to clarify the gap structure of 
iron-based superconductors^ 

In this Rapid Communication, we present the ther- 
mal conductivity measurements of optimally doped 
NaFeo.97 2 Coo.o28As and overdoped NaFeo.9 2 5Coo.o7sAs 
single crystals down to 50 mK. We find no evidences 
for gap nodes in both compounds, which supports pre- 
vious ARPES results in Ref. 23, and disagrees with 
the penetration depth experiments in Ref. 24. Fur- 
thermore, the relatively fast field dependence of kq/T 
in NaFeo.925Coo.o75As indicates the increase of the ra- 
tio between the magnitudes of those isotropic gaps, or 
the enhancement of gap anisotropy upon increasing dop- 
ing. This gap structure evolution is similar to that of 
Ba(Fei„ x Co x ) 2 As2 system. 

The NaFei-zCozAs single crystals were synthesized by 
flux method with NaAs as the fluxj^l NaAs was prepared 
by reacting the Na chunks and As powders in the evac- 
uated quartz tube at 200 °C for 10 hours. Then the 
powders of NaAs, Fe and Co were mixed together ac- 
cording to the ratio NaAs:Fe:Co = 4:l-x:cc. The mix- 
ture was placed into an alumina crucible and then sealed 
inside an iron crucible under Ar atmosphere. The cru- 
cible was put in Ar filled tube furnace and slowly heated 
up to 950 °C and kept for 10 hours, then slowly cooled 
down to 600 °C at a rate of 3 °C/h. The shining plate- 
like NaFei_ x Co x As single crystals with maximum size of 
5x5x0.2 mm 3 were obtained. The actual chemical com- 
position of the single crystals is determined by energy- 
dispersive X-ray spectroscopy. The cobalt contents are 
x = 0.028 and 0.075, respectively, for the samples with 
nominal x = 0.05 and 0.20. The dc magnetic susceptibil- 
ity was measured at H = 20 Oe, with zero-field cooled, 
using a SQUID (MPMS, Quantum Design). 

The single crystals were stored and shipped under the 
environment of inert gas. After exposed in air, the sam- 
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FIG. 1: (Color online), (a) Normalized dc magnetization 
of NaFeo.972Coo.o28As and NaFeo.925Coo.075 As single crys- 
tals used for transport study, (b) Magnetization of another 
NaFeo.925Coo.o75As single crystal from the same batch, which 
also has millimeter size in the afc-plane and ~100 /im along 
c-axis, after being exposed in air for 0.2, 2.2, and 12 hours, 
respectively. The magnetization barely changes during the 
initial 2.2 hours. 



FIG. 2: (Color online). (a) In-plane resistivity of 
NaFeo.972Coo.o2gAs and NaFeo.92sCoo.o75As single crystals. 
The solid lines are linear fit for NaFeo.972Coo.028 As be- 
tween 50 and 140 K, and a fit to p(T) = po + AT" for 
NaFeo.925Coo.o75As between 25 and 140 K, respectively, (b) 
The upper critical field Hca(T) of NaFeo.972Coo.o2sAs and 
NaFeo.925Coo.075 As, defined by p = 0. For comparison, the 
data of Ba(Feo.9Co .i)2As 2 with T c w 22 K is also plotted. 29 



pie was quickly cleaved to a rectangular shape with typ- 
ical dimensions of ^2.50x1.00 mm 2 in the a&-plane, and 
^100 /im along the c-axis. Then four silver wires were 
attached on the sample with silver paint, which were 
used for both thermal conductivity and resistivity mea- 
surements. The time of exposure in air is less than one 
hour, before the heat transport study. In-plane thermal 
conductivity was measured in a dilution fridge using a 
standard four-wire steady-state method with two Ru02 
chip thermometers, calibrated in situ against a reference 
RUO2 thermometer. After the thermal conductivity mea- 
surements, the sample was quickly (roughly 30 minutes 
exposure in air) switched to a 4 He cryostat for resistivity 
measurements. Finally, the contact resistances were ex- 
amined, with typical value of 100 mil at 2 K. Magnetic 
fields were applied along the c-axis. To ensure a homo- 
geneous field distribution in the samples, all fields were 
applied at temperature above T c . 

Figure 1(a) presents the normalized dc magnetiza- 
tion of NaFeo.972Coo.o28As and NaFeo.925Coo.o7sAs single 
crystals. The onset transition temperatures T c are 20 and 
11 K, respectively. For NaFeo.972Coo.o28As, the data were 
taken on a crystal right after being cleaved from the sam- 
ple used for transport study. For NaFeo.925Coo.o75As, the 
data were taken directly on the sample used for transport 
study, after all transport experiments have been done, 
totally about 2 hours exposure in air. In order to check 



how sensitive is the bulk superconductivity to the expo- 
sure time, we measured another NaFeo.925Coo.o7sAs sin- 
gle crystal from the same batch, which also has millimeter 
size in a6-plane and ^100 /im along c-axis. Fig. 1(b) 
shows its magnetization after being exposed in air for 
0.2, 2.2, and 12 hours, respectively. The magnetization 
barely changes during the initial 2.2 hours, therefore it 
is concluded that within 2.2 hours there is little effect of 
air on the bulk superconductivity of NaFei-^Co^As sin- 
gle crystals with this kind of size. We notice that the bulk 
superconductivity degrades faster in air for smaller and 
thinner single crystals, which is understandable. Since 
the contact resistances remain low at the end of trans- 
port measurements, our resistivity and thermal conduc- 
tivity data represent the intrinsic properties of the two 
NaFei^CozAs compounds. 

Figure 2(a) shows the in-plane resistivity p(T) of 
NaFe .972Co .o28As and NaFe .925Co .o75As single crys- 
tals. The T c defined by p = are 20.2 and 12.4 K, respec- 
tively. This T c of NaFeo.925Coo.o75As sample is about 1 
K higher than that determined from magnetization mea- 
surement, possibly due to slight inhomogeneity. We name 
these two samples as OP20K and OD11K. For OP20K, 
the resistivity data between 50 and 140 K can be fitted 
linearly, giving a residual resistivity po = 167.8 ± 0.1 
pClcm. For OD11K, the data between 25 and 140 K are 
fitted to p{T) = po + AT n , which gives p = 151.6 ± 0.1 
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FIG. 3: (Color online). Low-temperature in-plane thermal 
conductivity of NaFeo.972Coo.o2sAs and NaFeo.925Coo.075 As 
in zero and magnetic fields applied along the c-axis. The 
solid lines are fits to k/T = a + bT 01-1 . The dash lines are the 
normal-state Wiedemann- Franz law expectation Lo/po, with 



Lo the Lorenz number 2.45 x 
and 151.6 /if2cm, respectively. 



1CT 8 WfiK and p = 167.8 



/ificm and n — 1.76 ± 0.01. The non-Fermi-liquid linear 
behavior of p(T) near optimal doping, and the increase of 
power n in the overdoped regime have been observed in 
BaFe2(Asi_ a; P : i;)2, which was considered as the signature 
of a quantum critical points 

The resistivity of these two samples were also mea- 
sured in magnetic fields up to H = 12 T, in order 
to estimate the upper critical field H C 2- Fig. 2(b) 
plots the temperature dependence of H C 2 for OP20K and 
OD11K, defined by p = 0. For comparison, the data of 
Ba(Feo. 9 Coo.i)2As 2 with T c « 22 K is plotted^ From 
Fig. 2(b), we roughly estimate H C 2 ~ 36 and 15 T for 
OP20K and OD11K, respectively. To choose a slightly 
different H C 2 does not affect our discussion on the field 
dependence of kq/T below. 

Figure 3 shows the temperature dependence of in-plane 
thermal conductivity for OP20K and OD11K in zero and 
magnetic fields, plotted as k/T vs T. To get the residual 
linear term kq/T, we fit the curves to k/T — a+bT a ~ 1 , in 
which the two terms aT and bT a represent contributions 
from electrons and phonons, respectivel y 30 ' 31 The power 
a of the phonon term is typically between 2 and 3 i 30 i 31 
For OP20K in zero field, the fitting gives kq/T = a = 
0.004 ± 0.003 mW K~ 2 cm" 1 , with a = 2.26 ± 0.03. 
For OD11K in zero field, kq/T = -0.005 ± 0.008 mW 
K~ 2 cm^ 1 and a — 2.32 ± 0.06 are obtained. Previously 



a w 2.2 has been observed in BaFe2As2r^ and a 
found in BaFe1.9Nio.1As2 and KFe 2 As 2 i&± - 
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FIG. 4: (Color online). Normalized residual linear term 
ko/T of NaFeo.925Coo.o75As as a function of H/H C 2- For 
comparison, similar data are shown for the clean s-wave 
superconductor Nb,— the dirty s-wave superconducting al- 
loy InBi, 35 the multi-band s-wave superconductor NbSe2,— 
an overdoped d-wave cuprate superconductor Tl-2201, 33 
the optimally doped Ba(Feo.92eCoo.o74)2As2 and overdoped 
Ba(Feo.892Coo.ios)2As 2 .- 



Both kq/T of OP20K and OD11K in zero field are 
negligible, within our experimental error bar ± 0.005 
mW K~ 2 cm _1 4i In nodelss superconductors, all elec- 
trons become Cooper pairs as T — > 0, therefore there 
are no fermionic quasiparticles to conduct heat and zero 
kq/T is observed, as in the conventional s-wave super- 



conductor V3SL 3 - On the contrary, in a nodal supercon- 
ductor, the nodal quasiparticles will contribute a finite 
kq/T in zero field, as in the <i-wave cuprate superconduc- 
tor Tl 2 Ba2Cu0 6+d - (Tl-2201)^ 3 - Therefore the negligible 
kq/T in both OP20K and OD11K strongly suggest that 
their superconducting gaps are nodeless. These results 
are consistent with previous ARPES results in Ref. 23, 
but disagree with the London penetration depth experi- 
ments in Ref. 24. 

To gain further support for nodeless gap in OP20K 
and OD11K, we check the field dependence of their kq/T. 
For OP20K, the fit to the data in H = 3 T also gives a 
negligible k /T = 0.004 ± 0.011 mW K~ 2 cm" 1 . Fur- 
ther increasing field, the data in H = 6 and 9 T become 
noisy and the fits are not very good, but the tendency 
is that the kq/T does not increase noticeably up to H = 
9 T (« i/ c2 /4). This field effect on k /T of OP20K is 
again consistent with the ARPES experiments on opti- 
mally doped NaFei-zCo^As, which show isotropic gaps 
with Ai = 6.8 meV for the a hole band and A2 = 6.5 
meV for the ^(5) electron bands^ 3 - Since the magnitudes 
of these gaps are very close (Ai/A 2 = 1.05), the field 
effect on kq/T is just like that of a superconductor with 
single isotropic gap^i 3 - for examples, the clean s-wave su- 
perconductor Nb and the dirty s-wave superconducting 
alloy InBi, shown in Fig. 4j 34 i 35 

As seen in Fig. 3(b), the fits to the thermal con- 
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ductivity data of OD11K in fields look reliable, so all 
k /T are obtained up to H = 9 T. In Fig. 4, the 
normalized {kq/T)/(knq/T) of OD11K is plotted as a 
function of H/H C 2, with the normal-state Wiedemann- 
Franz law expectation knq/T = Lo/po = 0.162 mW K -2 
cm -1 and H C 2 = 15 T. Similar data of the multi-band s- 
wave superconductor NbSe 2 ^ Tl-2201^ the optimally 
doped Ba(Feo. 926Co . 074)2As2 (T c = 22.2 K) and over- 
doped Ba(Feo.892Coo.io8)2As2 (T c = 14.6 K)2 are also 
plotted for comparison. Clearly, the Kq(H)/T behav- 
ior of OD11K is very similar to that of optimally doped 
Ba(Feo.926Coo. 074)2 As2, and much slower than that of 
overdoped Ba(Fe . 892C00. 108)2 As 2 . 

For optimally doped BaFei.85Coo.isAs2, the gaps are 
isotropic, with Ai = 6.0 meV for the j3 hole pocket and 
A2 = 5.0 meV for the 7(c)) electron pockets^ The dif- 
ferent magnitudes of those gaps (A1/A2 = 1.32) may 
result in the relatively fast field dependence of kq/T 
shown in Fig. 4. For overdoped Ba(Feo. 892C00. 108)2 AS2, 
the k,q/T increases very rapidly even at low field, which 
has been interpreted as the manifestation of large ra- 
tio (> 3) between Ai and A2r^ or highly anisotropic 
gap<2 The similar kq(H)/T behaviors between OD11K 
and optimally doped Ba(Feo.926Coo.o74)2As2 suggest that 
the ratio A1/A2 has increased from 1.05 in OP20K to 
~ 1.3 in OD11K, or some gaps become anisotropic up 
to 30%. This needs to be checked by the ARPES ex- 
periments. From our data, the gap structure evolu- 
tion upon increasing doping is similar in NaFei_ x Co x As 



and Ba(Fei_ x Co x )2As2 systems. Note that the field de- 
pendence of kq/T for OD11K again disagrees with the 
claim of nodal gap in overdoped NaFei_ x Co^As,— since 
k (H)/T should display a much faster field dependence 
in that case. 

In summary, we have measured the thermal conductiv- 
ity of optimally doped NaFeo.972Coo.o2sAs and overdoped 
NaFeo.925Coo.o75As single crystals down to 50 mK. The 
absence of kq/T in zero field for both compounds gives 
strong evidence for nodeless gap. The field effect on 
Ko/T of NaFe .972Coo.o28As is consistent with multiple 
isotropic gaps with similar magnitudes, as demonstrated 
by ARPES experiments. The relatively faster field de- 
pendence of kq/T for overdoped NaFeo.925Coo.o7sAs sug- 
gests that the different between the magnitudes of those 
gaps increases, or some gaps become anisotropic upon 
increasing doping. Our results clearly disagree with the 
claim of nodal gap in NaFei-^Co^As based on London 
penetration depth experiments, therefore the issue of gap 
structure in this 111 system has been clarified. 
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